Absence of a true long-range orbital order in a two-leg Kondo ladder 
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We investigate, through the density-matrix renormalization group and the Lanczos technique, the 
possibility of a two-leg Kondo ladder present an incommensurate orbital order. Our results indicate a 
staggered short-range orbital order at half-filling. Away from half-filling our data are consistent with 
an incommensurate quasi-long-range orbital order. We also observed that an interaction between 
the localized spins enhances the rung-rung current correlations. 
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I. INTRODUCTION 

In 1985, it was observed that the heavy fermion su- 
perconductor URu 2 Si2 presents a second order phase 
transition at 17. 5K. 1 This phase transition is character- 
ized by sharp features in the specific healji and several 
others thermodynamic properties (see, e.g., Ref. and 
References therein). The large entropy loss associated 
in this phase transition is equivalent to an ordered mo- 
ment of about 0.5/xs. However, the size of staggered mo- 
ment measured by neutrons scattering measurements is 
m ~ The order parameter associated with this 

phase transition is, at the present moment, not estab- 
lished and it is challenging to discover the nature of the 
hidden order behind the transition. 

Many theoretical groups have proposed several kinds of 
hidden order But, until now, experiments were 
not able to establish which is the correct one. Certainly, 
also from the theoretical point of view, more studies are 
needed to clarify the correct order associated with this 
mysterious phase transition. In this front, we present 
here a numerical study of a microscopic model for the 
heavy fermion systems. 

In this work we focus on the order parameter proposed 
a few years ago by Chandra and collaborators.— They sug- 
gested the existence of a hidden incommensurate orbital 
order in the heavy fermion URU2S12 below the second 
order phase transition. The orbital order phase is associ- 
ated with currents circulating around the plaquettes, as 
illustrated in Fig. 1. In the case of URu2Si2, this cur- 
rents produce a very week orbital moment 0.02/is that 
explains the large entropy Iossj^ 

Very recently, neutron scattering measurements were 
unable to detect the orbital order in the heavy fermion 
URu2Si2- 8 Although the orbital order was not detected 
it is not possible yet to discard it as the hidden order due 
to the resolution limitation of the experiments performed. 
Note that the orbital order is expected to be 50 times 
smaller than the spin order £ 

Our goal in this work is to investigate the existence of 
an incommensurate orbital order in the Kondo Lattice 
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Figure 1: (Color online) A schematic representation of the 
two-leg Kondo ladder. It also shows the circulating currents 
around the plaquettes (in this example a staggered one). 



model (KLM). This model is the simplest one believed 
to present the physics of heavy fermions materials^ (see 
next section). Our approach will be numerical, through 
the density-matrix renormalization group (DMRG)^ and 
the Lanczos technique. 12 These techniques are non- 
perturbative, however limited by the system size. For 
this reason, we consider the two-leg Kondo ladder (2- 
LKL), which is the simplest geometry able to present an 
orbital order. 

The orbital order, also called flux or orbital current 
phase, has already been discussed in the context of the 
high temperature superconductors. The standard two-leg 
t-J ladders model present a short-range orbital order^ 
while an extended version has long-range orbital order 
for some parameters^ A recent detailed discussion of 
the orbital order in the context of a Hubbard model can 
be found in Ref. [FoI . 

We close this section mentioning that a model very sim- 
ilar to the KLM was used to describe the magnetism of 
URu2Si2- Sikkema and collaborators^, through a mean 
field calculation, showed that the Ising-Kondo lattice 
model with transverse field presents a weak ordered mo- 
ment, similar to the one observed in experiments. How- 
ever, the Ising-KLM model was not able to reproduce the 
large specific heat jump. 
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II. MODEL 

In order to investigate the heavy fermion systems the 
minimum ingredients that a microscopic model must con- 
sider are two types of electrons, the conduction electrons 
in the s-, p-, or d-, orbital as well the electrons in the 
inner f-orbitals.— In the literature there are two well 
known standard models that consider these two kind of 
electrons, the periodic Anderson model (PAM) and the 
KLM. 17 In an appropriate parameter regime (mainly (i) 
the mobility of the / electrons is very small, which is 
relevant for the heavy fermion system and (ii) that the 
Coulomb interaction of the electrons in the / orbitals 
is very large) Schrieffer and Wolf f 17 ' 18 showed that the 
KLM can be derived from the PAM. We consider in this 
work the KLM which has less degrees of freedom per unit 
cell than the PAM and it is easier to explore numerically. 

The KLM incorporate an interaction between the local- 
ized spins and the conduction electrons via exchange in- 
teraction J. To attack this model in two or three dimen- 
sion by unbiased non-perturbative numerical approaches 
is an impossible task at the present moment. However, 
it is possible to consider quasi-one-dimcnsional systems 
such as the N-leg ladders model. 

We consider the 2-LKL with 2xL sites defined by 

H KM = - E (4, CT S>+ H - c -) + J E S i' s ^ 

<i,j>,a j 

+ Jam 22 Si • Sj, (1) 

<i,j> 

where Cj a annihilates a conduction electron in site j with 
spin projection a, Sj is a localized spin ^ operator, Sj = 

\ Ea/3 c ] a <T a/3 c j B ' ls the conduction electron spin density 
operator and <r a g are Pauli matrices. Here < ij > denote 
nearest-neighbor sites, J > (when the KLM is deduced 
from the PAM obtain J > 0i£) is the Kondo coupling 
constant between the conduction electrons and the local 
moments and the hopping amplitude was set to unity to 
fix the energy scale. 

We also consider an interaction between the localized 
spins Jah, we choose Jam > since antiferromagnetims 
had been observed in URu2Si2^ The same model above 
also represents the manganites when J < 0i2°- In this lat- 
ter case, the interaction between the localized spins seems 
to be important to stabilize some phases^ - This is the 
motivation to also consider this interaction. Note that 
several others terms in the Hamiltonian could also be in- 
cluded, like the Coulomb interaction of the electrons in 
the conduction band, extra electrons hopping, etc. How- 
ever, at the present moment, there are no evidences in- 
dicating that such extra terms are relevant to the low 
energy physics of the heavy fermion systems. Up to now, 
it is well established that J is essential to describe the 
magnetism observed in the heavy fermion systems. At 
small values of J, an antiferromagnetic long-range or- 



der (LRO) is expected due the Ruderman-Kittel-Kasuya- 
Yosida interaction, whereas for large J a paramagnetic 
phase emerges. Doniach^i was the first to point out the 
existence of a quantum critical point (QCP) due the com- 
petition between these two phases. 

Unlike other models, such as the t — J model, much less 
is known about the Kondo lattice model. Even in the one 
dimension version, where the ground state of the Kondo 
chain is quite well known 1 ^ (see also Ref. l22l ). New 
phases have been reported recently, such as a new ferro- 
magnetic phased inserted into the paramagnetic phase 
as well as a dimerized phase at quarter-filling. 24 The lat- 
ter has been questioned recently by Hotta and Shibata. 25 
Those authors claim that the dimerized phase is an arti- 
fact of the open boundary conditions. Indeed, the bound- 
ary condition is very important, as well as the number of 
sites considered. In Ref. H^the authors observed, mainly, 
that with an odd number the sites the dimer state does 
not exist. The parity of the number of sites is thus very 
relevant and an odd number destroys the dimerization. 26 

In quasi-one-dimensional systems, such as the N-leg 
ladders, very few non-perturbative studies have been re- 
ported. Recently, quantum Monte Carlo 2 ^ (QMC) and 
DMRG 28 calculations of the half-filled Kondo lattice 
model in small clusters found the existence of a quantum 
critical point (QCP) at J ~ 1.45, in agreement with pre- 
vious approximated approache s 29 ' 30 (see also Ref. [3lh - 
Note that the QMC calculations were feasible only at 
half- filling, where the famous sign problem is absent. 
Moreover, the DMRG results of the N-LKL at half-filling 
show that the spin and charge gaps are nonzero for any 
number of legs and Kondo coupling J. These results 
are quite different from the well known N-leg Heisenberg 
ladders were the spin gap is zero for an even number of 
legs i^ 2 - 

The phase diagram of the 2-LKL has also been ex- 
plored numerically— In this CcLS6, cL ferromagnetic phase 
was observed only for small densities, very distinctively 
from the phase diagram of the ID Kondo lattice chain, 
where the ferromagnetism is present at all electronic den- 
sities for large values of J. However, it is similar to 
the mean field phase diagram of the 3D Kondo lattice 
models In this sense, the 2-LKL presents a better sig- 
nature of the phases appearing in real systems than its 
one-dimensional version. Interesting that it was also ob- 
served dimerization in the 2-LKL^ 3 - at conduction elec- 
tron densities n = 1/4 and n = 1/2. As in the one- 
dimension version, the RKKY interaction explains these 
unusual spin structures. In fact, in some real heavy 
fermion systems some unusual spin order structures have 
indeed been observed^ 

Here, we consider electronic densities n larger than 0.4, 
where a paramagnetic phase have been observed^ In 
particular, we focus on the electronic densities n = 1 and 
n = 0.8. We choose these densities since the magnitude 
of the rung-rung current correlation is bigger for larger 
electronic densities. We investigate the model with the 
DMRG technique under open boundary conditions and 
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use the finite-size algorithm for sizes up to 2 x L = 120, 
keeping up to m = 1600 states per block in the final 
sweep. The discarded weight was typically about 10~ 5 — 
10~ 7 in the final sweep. We also cross-checked our results 
with Lanczos technique for small systems. 

III. RESULTS 

Before presenting our results, we briefly discuss the 
order parameter associated with a circulating current 
phase. Such a phase breaks rotational, translational as 
well as time reversal symmetries. The appropriated or- 
der parameter to detect this phase is the current between 

two nearest-neighbour sites, i. e., (^Jl,jj where the cur- 
rent operator between two nearest-neighbours i and j is 
given by 

Strictly speaking, a spontaneous symmetry breaking 
only appears in the thermodynamic limit. Only in this 

limit m the ordered phase. The signature 

of a spontaneous symmetry breaking appears in the two 
point correlation function of the operator that measures 
the symmetry. We utilize this fact to infer about the or- 
bital order. If a continuous symmetry is broken, no long- 
range order exist at finite temperature in one and two 
dimensions, as stated by the Mermin-Wagner-Hohenberg 
theorem.— At zero temperature a true long-range order 
is still possible in two dimensions, while in one dimen- 
sion only a quasi-long-range order can manifest, i. c., the 
two point correlation function decay algebraic. However, 
if the translational symmetry (a discrete symmetry) is 
broken, even in one dimension a true long-range order 
may exist (a famous example is the dimerized phase of 
the Majumdar- Ghosh mode l 37 i 38 ' 39 ). Since the transla- 
tional symmetry is broken in the orbital phase, a true 
long-range order may occurs in the ground state of the 
2-LKL. 

In order to observe any trace of orbital order in the 
ground state wave function of the 2-LKL, we measure 
the rung-rung current correlations defined as 

C(l,k) = (j(l)J(k)), 

where J(l) is the rung current operator for the Zth rung 
given by 

J{l)=i X)( c k<7 C !i,<7 ~ 4 1(T ci 2}tT ), (2) 

and c;a,o- annihilates a conduction electron on rung I and 
leg A = 1, 2 with spins a =||. Since we work with open 



boundary conditions it is convenient to define an aver- 
aged rung-rung current correlation in order to minimize 
boundary effects. We have defined the averaged rung- 
rung current correlation as 

C ^ = ji E (-mm), (3) 

\i-k\=l 

where M is the number of site pairs (i, k) satisfying / = 
\i — fc|. Typically, M in our calculation vary from 3 to 
10. 

There is a true long-range orbital order if 
limi^ 00 C{l) 7^ 0. Through this criterion, we can 
infer about the existence of the orbital order by measur- 
ing C(l) at large distances. If C(l) has an exponential 
decay, the linear-log plot shows a linear decay. On the 
other hand, if C(l) has an power law decay, the log-log 
plot present a linear decay. 

We also measure the cosine transform of C(l), i. e., 

L 

N(q) = ^2C(l)cos(lq), 
i=i 

in order to infer about periodicity of the oscillatory part 
ofC(Z). 

In the next two subsections we investigate these corre- 
lations for the 2-LKL at half-filling and close to half-filing, 
respectively. We did not find any evidence of long-range 
orbital order in the ground state of the 2-LKL. Our re- 
sults support that the rung-rung current correlation has 
an exponential decay at half-filling. Close to half-filling 
our results indicate an incommensurate quasi-long range 
orbital order. 



A. Half-filling 

We start presenting some results for the conduction 
density n = 1. We observed, in this case, that the aver- 
aged rung-rung current correlation behaves as 

C(0=a o (-l) I expH/O, (4) 

for all values of J and J ah explored in this work. In Fig. 
2(a), we present a typical example of the magnitude of 
C(l) at half-filling for a system size L — 30. As we see, 
our results indicate strongly that C(l) has an exponential 
decay due to the linear decay in the linear- log plot. The 
inset in Fig. 2(a) also shows that C(l) is staggered. The 
solid line in Fig. 2(a) correspond to a fit of Eq. 2 with 
ao = 0.16 and a decay length £ = 1.434^ We performed 
a least-squares fitting, resulting in a root mean square 
(RMS) of 0.0018 and a correlation coefficient of 0.996. 
We found that C(l) has a very small dependence on the 
number m of states retained in the truncation process for 
J > 0.8, as can be observed in Fig. 2(a). For J < 0.8 is 
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Figure 2: (Color online) (a) The linear-log plot of |C(Z)| for 
two distinct value of m with L = 30 at half-filling. The solid 
line in Fig. 2(a) correspond to a fit of Eq. 2 with £ = 1.43 
and ao — 0.06, the RMS per cent error is 0.18. Inset: C(l) vs 
distance with m = 1000. Only few sites are presented. The 
couplings are J = 0.8 and J ah = 0. (b) The cosine transform 
N(q) of C(l) presented in Fig. 1(a) with m = 1000. (c) 
The linear-log plot of \C(l)\ for two distinct size, both with 
m — 1000. The couplings are J = 0.35 and J ah = 0. 



very hard to get accurate results, however even for small 
J we believe to have captured the correct qualitative be- 
haviour. Nevertheless, we present most of our results for 
J > 0.8, where the results are more accurate. 



Figure 3: (Color online) The linear-log plot of \C(l)\ for a 
set of representative values of J and J ah for L = 30 at half- 
filling, (a) \C(l)\ for J ah = and J = 0.35, 0.8 and 1.8. (b) 
\C(l)\ for J = 0.35 and some values of Jah- (c) Same as (b) 
but for J = 1.8. 



The signature of the sign alternation is observed 
through the cosine transform of C(l). In Fig. 2(b), we 
show the cosine transform of C(l) present in Fig. 2(a) 
with m = 1000. Clearly, we observed a peak at q = it 
due the sign alternation of C(l). Note that the finite-size 
effects are small, as can be seen in Fig. 2(c). For this 
reason, we restrict most of our calculus to system size 
2x30 in order to save computational time. 
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Figure 4: (Color online) The peak intensity of N(q) at q — tv 
for J — 0.8 and J = 1.8 as a function of J ah for a system 
size L = 30 at half-filling. 



Our results indicate that for small J, where the RKKY 
is expected to be dominant, the rung-rung current corre- 
lations has a bigger correlation length as we see in Fig. 
3(a). On the other hand, for large J, which favors for- 
mation of singles, the correlation length is smaller. This 
result is expected, since for J — * oo the rung-rung current 
correlations must go to zero. 

The Hamiltonian Eq. (1) with J ah = does not lead 
to a long-range orbital order at half-filling, as we have 
observed. Since J ah seems to be important to stabilise 
some phases for J < 0^2 ; it may be possible that it also 
stabilises the orbital phase for J > 0. For these reason, 
we also investigate the effect of J ah in the ground state 
of the 2-LKL. As we see in Fig. 3(b), for small values of 
J, J ah does not affect significantly C(l). On the other 
hand, for larger J as shown in Fig. 3(c), J ah clearly en- 
hances the length correlation. Although J ah enhanced 
C(l), at half-filling only short-range orbital order is ob- 
served for several parameters investigated. 

At half-filling, for all parameters studied, N(q) always 
presents a peak at q = n. In Fig. 4, we present this 
peak intensity for J = 0.8 and J = 1.8 as function of 
J ah- As we see, the peak intensity increases with J ah 
and saturates for large J ah around ~ 0.45. 

Our main conclusion, for the half-filling case, is absence 
of long-range orbital order. Note that it may be possible 
that the inclusion of the Coulomb interaction between 
the electrons in the conduction band leads the system 
into a phase with long-range orbital order, as occurs in 
an extend t — J model.— This is under investigation in 
the present moment by one of the authors. 



B. Close to half-filling 

Away from half- filling the DMRG calculation of C(l) 
is less stable, for this reason we consider system sizes 
smaller than 2x40 and keeping up to m = 1600 states in 
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Figure 5: (Color online) The log-log plot of \C(l)\ for the 
conduction density n = 0.8. (a) \C(l)\ vs distance for J = 0.8, 
J ah = 1.0, and L = 30 for two distinct value of m. The 
dashed line is the fit using Eq. (1) with ai = 2, 012 = 3, 
do — 0.16 and ai = —0.33, the RMS per cent error is 6 with 
a correlation coefficient of 0.998. The solid line is guide by 
eyes, (b) |C(£)| for the same parameters of Fig. 5(a) and 
two distinct sizes, both cases with m = 1600. Inset: C(l) vs 
distance with m = 1600. Only few sites are presented. 



the truncation process. Although we obtained results for 
a few densities away from half-filling, we focus on density 
n = 0.8 where the magnitude of C(l) is larger. For small 
densities is very hard to get accurate results since the 
current intensity is very small. In Fig. 5(a), we present 
the log-log plot of \C(l) \ at conduction density n = 0.8 for 
a system size 2 x 30 with J = 0.8 and J ah = 1.0 for two 
different values of m. Since in the log-log plot we obtain 
a linear decay (see the solid line in this figure) C(l) must 
have a power law decay. If we use a linear-log plot our 
data does not have a linear decay. As can be seen from 
that Figure, it is very hard to get good accuracy even 
working with m = 1600 states. Although we were not 
able to obtain the current-current correlations at large 
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kept in the truncation process since we also obtained the 
same effect for small clusters with exact diagonalization. 

In Fig. 6 we present N(q) for a representative set of 
parameters at conduction density n = 0.8. As shown 
in that Figure, there is no peak at q — tt, signaling an 
absence of staggered rung-rung current correlations. For 
the conduction density n = 0.8 we observed a cusp at 
q = nn. These results indicate that close to half-filling 
the 2-LKL presents an incommensurate quasi-long-range 
orbital order. 



IV. CONCLUSION 



Figure 6: The cosine transform N(q) of C(l) presented in 
Figure 5(a) with m = 1600. 



distances with a high accuracy, we believe to have cap- 
tured the correct behavior, i.e., a power law decay. The 
large oscillations appearing in those Figures are due to 
fact that some values of C(l) are very close to zero. 

Since our data of C(l) in the log-log plot strongly sug- 
gest a power law decay close to half-filling (note that for 
the half-filling case the decay is exponential) we tried to 
fit C(l) with the function 



„ ... cos(mrl) cos(2nirl) 
C fit (l)=a — ±—L+ ai - 



la 2 



(5) 



where n = 0.8 is the density. The dashed curve in Fig. 
5(a) corresponds to a fitting of our data with m = 1600. 
We were not able to reproduce precisely C(l), however 
the general behavior is quite well described. 

Note also that finite-size effect are larger away from 
half-filing, as we can see by comparing the Figs. 5(b) and 
2(c). It is important to mention that we observed, away 
from half- filling and in very few distances I, that the sign 
of the averaged correlation C(l=|j-k|) does not has the 
same sign of C(j, k) for some pairs of (j,k) satisfying l=|j- 
k|. This does not seem to be due to the number of states 



In this paper, we have investigated the possibility of a 
two-leg Kondo ladder present an orbital order. In partic- 
ular, we focus on the densities n = 1 and n — 0.8. For the 
several couplings investigated we did not find any trace of 
a true long-range orbital order, which would be relevant 
to explain the large entropy loss observed in the second 
order phase transition of URu2Si2 ■ Our data indicate 
that the half-filling case presents a staggered short-range 
orbital order, while close to half-filling our results are con- 
sistent with an incommensurate quasi long-range orbital 
order. Although we did not find evidence of a long-range 
orbital order in the ground state of the two-leg Kondo 
ladder, we can not yet completely discard this possibil- 
ity. It may occur that an extended version of the Kondo 
lattice model presents the long-range orbital order. So, 
we may conclude that either the orbital phase does not 
exist and is not the origin of the mysterious phase tran- 
sition observed in the the heavy fermion U RU2S12 or the 
standard Kondo lattice Model is not able to reproduce 
the correct order observed in the experiments. 
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